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摘 要 : 两 相 流 系 统 中 ， 流 型 会 直接 影响 系统 的 摩擦 阻力 和 传 热 等 特性 ， 而 流 型 的 转变 过 程 会 受到 通道 的 形状 和 尺寸 
的 影响 。 常 温 常 压 下 ， 以 空气 和 水 为 介质 ， 对 八 种 不 同 尺 寸 矩形 通道 中 竖 直 向 上 泡 状 流向 弹 状 流转 变 过 程 进行 了 下 
究 。 结 果 表 明 ， 拢 形 通道 间隙 大 于 2 mm 时 ， 随 着 通道 水 力 直 径 的 增 大 ， 泡 状 流向 弹 状 流转 换 边界 左 移 。 在 窗 通 道 范 
围 内 ， 可 采用 Ishii (1977) 提出 的 基于 漂移 流 模型 的 判定 准则 计算 流 型 转换 边界 ， 在 常规 矩形 通道 中 ， 转 换 准 则 中 的 
临界 空 泡 份额 与 气泡 初始 尺寸 有 关 ， 选 用 赵 建 福 提出 的 方法 计算 临界 空 泡 份额 ， 转 换 边 界 预测 结果 与 实验 结果 符合 程 
度 较 好 ， 在 利用 漂移 流 方法 计算 截面 含 气 率 时 ， 对 于 窗 矩 形 通 道 ， 分 布 系数 可 以 用 Ishii (1977) 提出 的 方法 计算 ， 而 
分 布 系数 为 1.2。 

关键 词 : 流 型 转换 ， 泡 状 流 ， 弹 状 流 ， 空 泡 份 额 ， 初 始 尺寸 
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Criterion for Predicting Transition Boundary from Bubbly Flow to Slug 


Flow in Rectangular Channels 
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Abstract: Flow pattern transition always affect the frictional resistance and heat transfer in two-phase 
systems, whereas it itself is also affected by the shape or size of the channels. The transition from bubbly 
flow to slug flow in eight different rectangular channels was investigated under normal condition. The 
results shows that the transition boundary from bubbly flow to slug flow moves to left as the diameter 
increasing when the channel gap is smaller than 2 mm. The transition criteria can be predicted from the 
model from Ishii (1977) in narrow rectangular channels. While for common channels, the critical void 
fraction in the transition criteria is relevant to the size of bubbles in entrance, for which calculation method 
of Jianfu Zhao (2005) is proposed, leading to the predicted transition criteria with a good agreement 
against the experimental result. When using of drift-flux model, the distribution parameter C, in channels is 
calculated with different forms, and the result shows that correction from Ishii (1977) can be used in 
narrow size of rectangular channels, while Co remains to be 1.2 in common channels. 
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两 相 流 系统 中 ， 流 型 的 转换 会 直接 影响 流动 阻 泡 状 流 、 
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Sadatomi 等 [对 水 力 直径 超过 10 mm 的 竖 


直上 升 矩 形 通道 内 流 型 进行 研究 ， 发 现 通 道中 存在 


弹 状 流 和 环 状 流 ， 当 直径 超过 10 mm 时 ， 


:特性 的 机 理 ， 基 于 流 型 的 经 验 公式 或 半 经 管道 的 几何 形状 对 流 型 的 影响 很 小 ，Mishima 等 外 
伟 公 式 也 将 随 之 发 生变 化 ， 因 此 研究 各 种 通道 中 的 通过 X 光照 相 方式 对 罕 矩 形 通道 中 的 流 型 进行 了 研 


Eh pire 都 要 以 精确 地 流 型 判别 为 基 完 ， 在 间隙 为 1.0 mm 的 鹤 矩 形 通 道中 搅 浑 流 没有 


础 。 对 于 和 矩形 通道 中 流 型 的 研究 主要 在 国外 ， 见 文 出 现 。Wilmarth 和 Ishii 等 人 BJ 对 间隙 为 1 mm 和 
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2mm HK FERE, JF 5 IshiibURI Taitellg 等 关 
于 泡 状 流向 弹 状 流 过 渡 的 模型 进行 验证 ， 并 提出 新 
的 用 于 计算 的 分 布 系数 Co。Hibiki 554188 3:158] E] [ft 
为 0.3 一 17 mm 的 矩形 通道 流 型 进行 研究 ， 得 到 了 
新 的 关于 流 型 转换 的 准则 。 

对 于 和 矩形 通道 内 流 型 转变 准则 的 研究 己 开 展 多 
年 ， 但 在 该 领域 得 到 的 结论 ， 尤 其 是 泡 状 流向 弹 状 
流转 变 的 准则 有 些 文献 给 出 的 结论 是 不 一 致 的 ， 本 
实验 选取 不 同 几何 尺寸 的 矩形 通道 ， 对 于 竖 直 向 上 
泡 状 流向 弹 状 流转 换 边 界 进行 研究 。 

1 实验 装置 

实验 在 常温 、 常 压 下 进行 ， 图 1 中 的 实 线 和 虚 
线 分 别 表示 水 回路 和 气 回 路 。 实 验 段 分 别 为 截面 为 
43x1.6mm? 〈 通 道 一 ) 、43x3mm? (通道 二 ) 和 
43xl0mm? GHiÉ —) 三 种 矩形 通道 ， 水 力 直径 分 别 
Jj 3mm. 5.6mm 和 16.2mm， 长 度 均 为 2000 mm, 
基于 KandlikarD] 的 定义 ， 通 道 一 为 罕 矩 形 通道 ， 通 
道 二 、 三 为 常规 矩形 通道 。 去 离子 水 由 离心 泵 提供 
驱动 力 ， 经 过 滤器 、 调 节 阀 和 混合 腔 进入 实验 段 返 
回 水 箱 形 成 循环 。 空 气 经 气泵 压缩 后 存在 储 气 钳 
中 ， 使 用 时 经 过 减 压 阀 和 止 逆 阀 进入 混合 腔 、 实 验 
段 后 由 实验 段 上 方 汽水 分 离 装 置 进入 大 气 。 混 合 腔 
中 均 布 若干 直径 为 1 mm 的 毛细 管 ， 确 保 实 验 段 入 
口 处 气泡 大 小 接近 ， 发 泡 均 匀 。 水 流量 和 气流 量 分 
别 采 用 不 同 质量 流量 计 测 得 ， 精 度 分 别 为 0.1 级 和 
0.2 级 。 距 实验 段 入 口 500 mm 及 1500 mm 处 安装 
有 精度 为 0.2 级 的 压力 传感器 。 正 对 实验 段 两 测 压 
点 中 心 位 置 固定 安装 高 速 摄影 仪 ， 用 以 拍摄 气泡 行 
为 。 
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图 1 实验 系统 图 


Fig.1 Schematic diagram of the experimental system 
2 泡 状 流 、 弹 状 流 流 型 的 判别 
本 实验 采用 高 速 摄影 方式 对 气 水 形态 


进行 观 
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测 ， 弥 补 人 眼 对 于 高 速 流体 辨识 不 清 的 缺点 。 文 献 
[1-4] 中 对 于 常规 矩形 通道 和 窜 矩 形 通道 中 的 泡 状 流 
和 弹 状 流 的 描述 是 不 同 的 ， 图 2 中 给 出 了 实验 中 不 
同和 矩形 通道 中 泡 状 流 和 弹 状 流 的 示意 图 。 

在 常规 通道 中 ， 泡 状 流 内 气泡 是 以 彼此 分 离 的 
小 泡 形式 分 散在 液体 之 中 ， 在 较 小 间隙 的 通道 中 ， 
气泡 多 为 圆柱 形 ， 较 大 间隙 的 通道 则 为 球状 ， 弹 状 
流 的 气 弹 多 数 有 圆润 的 头 部 和 局 平 的 主体 ， 气 弹 和 
液 相交 叉 分 布 ， 液 相 中 常 含有 大 量 弥散 的 小 气泡 。 
在 鹤 矩 形 通道 中 ， 由 于 通道 间隙 的 限制 ， 泡 状 流 即 
使 在 气相 速度 较 低 时 也 难以 发 现 ， 气 泡 提 前 聚合 形 
成 帽 状 流 的 流 型 ， 气 泡 头 部 呈 半 圆 型 ， 下 部 近似 为 
水 平 状 ， 弹 状 流 与 常规 矩形 通道 近似 相同 ， 但 窗 通 
道 弹 状 流 的 液 相 部 分 内 一 般 仅 含有 少量 的 气泡 ， 有 
时 气 弹 中 会 有 液 滴 出 现 。 


泡 状 流 AR 弹 状 流 


泡 状 流 RA 泡 状 流 
通道 一 通道 二 通道 三 
CD=3 mm) (D=5.6 mm) CD=16.2 mm) 


图 2 不 同 矩 形 通道 两 种 流 型 示意 图 


Fig.2 Two flow patterns in different kinds of rectangular 


channels 
3 泡 状 流 一 弹 状 流转 换 边 界 判定 方法 
3. 1 泡 状 流 一 弹 状 流转 换 边 界 对 比 


选取 从 文献 2 和 文献 11 得 到 的 和 珑 形 通道 中 泡 
状 流向 弹 状 流转 变 边界 ， 并 与 实验 中 得 到 的 结果 进 
行 比较 ， 如 图 3 所 示 《〈 图 中 s 为 通道 间 际 ，w 为 宽 
UKE, D 为 当量 直径 ， 尺 寸 单 位 均 为 mm)。 对 比 
发 现 矩形 通道 中 泡 状 流向 弹 状 流转 变 主 要 与 通道 间 
险 有 关 ， 当 通道 间隙 小 于 2 ”mm 时 ， 间 隐 为 0.3、 
1.0 和 1.6 mm 的 转换 边界 几乎 重合 ， 当 通道 间隙 大 
T 2mm 时 ， 随 着 通道 水 力 直径 的 增加 ， 转 换 边 界 
左 移 ， 由 于 间隙 变 小 造成 的 束缚 减弱 ， 壁 面 的 剪 切 
作用 也 相应 变 小 ， 气 泡 间 更 容易 实现 磁 撞 和 到 合 ， 


“ WORD 批量 转 PDF 工 具 - 未 注册 
注册 码 购买 耿 系 QQ:3049816538 


NN ` 日 二 和 
弹 状 流出 现时 间 相 应 提前 。 
10 
E 
dn 
P ， 
T 
ti 
ELUL 
K 
EX —TM— 5s—0.3mm,w-12mm,D-0.59mm 
i— s-0.6mm,w-12mm,D-71.14mm 
04 A s-].0mm,w-12mm,D-1.85mm 
4— s-].6mm,w-43mm,D-73.09mm 
—Ó— s-2.4mm,w-40mm,D-74.53nm 
—O-— s-3.0mm,w-43mm,D-75.61mm 
3X— s-5.0mm,w-40mm,/-8.89mm 
X— s-]0mm,w-43mm,D-16.2mm 
0.01 


0.01 0.1 


气相 折算 速度 (m/s) k 
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Fig.3 Comparison of different corrections for distribution 


parameter in rectangular channels 
3. 2 761A i S8 JA c8 38020 02] x EE DU HAY 
3. 2. 1 泡 状 流 一 弹 状 流转 换 边界 判定 准则 
许多 学 者 都 采用 漂移 流 的 公式 作为 泡 状 流向 弹 
状 流转 变 的 判别 式 ，Jones 和 Zuberls] (1975) 提出 
基于 矩形 通道 的 转换 边界 判定 准则 ， 形 式 为 
] ] 


Ją =- 0.23*0.13s / w J^pgw! p, 


7C C, 


(1) 

Ishii (1977〉 得 到 了 和 矩形 通道 泡 状 流 z% 

状 流 转换 准则 ， 指 出 气泡 直径 大 小 Dp 影响 了 临界 空 
泡 份额 的 值 : 


jC o ja IA u 
aC, C 


0 A 


(2) 


其 中 ， a —0.224$ « D , 


7 40.15% s 3D, 


0 = 
b 


a — 0.32437 3D, 
Taitel (1980) [I Sadatomi (1982) 上 相继 提 
出 了 适用 范围 较为 宽泛 的 转换 关系 式 ， 形 式 上 分 别 
为 


t=0.56j.—1.15(wgAp/ p? "^ (Taite) (3) 
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03C, 
可 以 看 出 ， 基 于 漂移 流 模型 提出 的 转换 准则 大 

部 分 含有 漂移 流 分 布 系数 Co， 因此 需要 得 到 不 同 矩 

形 通 道内 分 布 系数 的 值 。 

3. 2. 2 矩形 通道 内 分 布 系数 的 计算 
漂移 流 的 形式 一 般 为 ; 


(Je) (2) = 6 U)* Y. (5) 


其 中 <je> 为 平均 气相 折算 速度 ，< a> 为 平均 空 
泡 份 额 ，< 户 为 折算 速度 ， 到 ;为 平均 漂移 速度 。 
Zuber 和 Findley[9] 较 早 提出 了 关于 分 布 系数 的 计 


算 形式 Cu —1.2-0.24 p, / p, ， 随 后 Jones 和 


gj; JD, (Sadatomi) (4) 


Zuberi*!, Ishii14) JI psp F EE RUOBDÉ PEL Y 237 
系数 的 计算 形式 : Cr1.2 (Jones and Zuber) 和 


C, —1.35-0.35,] p, / p, (Ishii) 


实验 中 得 到 的 实验 段 的 总 压 降 AP 由 重 位 压 
降 、 摩 擦 压 降 和 加 速 压 降 组 成 。 由 于 通道 是 等 截面 
的 ， 在 绝热 条 件 下 ， 加 速 压 降 可 以 忽略 ， 摩 擦 压 降 
APr 经 由 Chisholm CL 方法 计算 得 到 ， 则 重 位 压 降 
为 


A © 
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即 : 


AP =| p1- (a) +p, (a) |gL (7) 


EH, <IRE ARIES 2 n UH LA 
两 测 压 点 间距 ， 于 是 利用 压 降 法 计算 出 的 平均 截面 
合 气 率 为 


AP- AP 
CT IE 
(a= & (8) 
Pi 7 Pe 


实验 得 到 的 结果 如 图 4 所 示 ， 随 着 矩形 通道 间 
隙 的 减少 ， 漂 移 流 模型 的 分 布 系数 有 逐渐 增加 的 趋 
势 ， 在 罕 通 道 范 围 ， 分 布 系数 可 以 选用 Ishii 进行 计 
算 ， 当 通道 间隙 大 于 5 mm 时 ， 分 布 系数 可 以 近似 
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Fig. 4 Comparison of distribution parameter in rectangular 


channels 
3. 2. 3 泡 状 流 弹 状 流转 换 准 则 评价 
根据 得 到 分 布 系数 的 计算 方法 ， 对 四 种 转换 准 
则 进行 计算 ， 误 差分 析 结 果 如 表 1 Hp. EPE 
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好 ， 这 主要 原因 是 在 罕 矩 形 通道 的 实验 研究 中 ， 气 
泡 初始 尺寸 基本 一 致 ， 对 于 流 型 转变 容易 形成 共 
识 ;， 而 在 较 大 间隙 的 通道 ， 由 于 气 液 混合 方式 不 一 
致 ， 造 成 气泡 初始 尺寸 相差 较 大 ， 模 型 一 认为 临界 
空 泡 份额 仅 在 0.2-0.3 变化 ， 因 此 有 必要 在 较 大 间 际 
的 矩形 通道 范围 对 Ishii 模型 的 临界 空 泡 份额 计算 式 
进行 修正 。 
文献 [12] 通 过 对 气泡 碰撞 对 于 泡 一 一 弹 状 流转 
换 进 行 数值 计算 ， 得 到 了 一 条 通用 的 、 关 于 气泡 初 
始 尺寸 影响 气泡 碰撞 率 的 通用 曲线 ， 研 究 表明 当 7 


C 2-DyD) 在 0.04-0.4 范围 内 ， 临 界 空 泡 份额 与 气 
泡 初 始 尺 寸 的 关系 是 
0 —1.18[(0.664 + 0.8767) ' - 1] (9) 


将 通过 高 速 摄影 捕捉 到 通道 二 与 通道 三 的 气泡 
初始 尺寸 应 用 到 公式 C 中， 并 用 Ishi 的 转换 准 
则 关系 式 的 公式 进行 计算 ， 得 到 新 的 泡 一 一 弹 状 流 
转换 边界 ， 与 实验 值 的 对 比 情况 如 图 5 所 示 ， 实 验 


形 通 道内 泡 状 流向 弹 状 流转 换 准 则 的 计算 中 ， 对 尺 
寸 较 小 的 通道 Ishii 模型 的 平均 误差 比较 小 ， 因 此 
Ishii 模型 可 以 用 于 窄 抢 形 通道 ， 而 对 于 较 大 尺寸 的 
和 矩形 通道 ， 四 种 模型 的 误差 都 较 大 ， 适 用 性 不 高 。 
Ishii 模型 通过 比较 气泡 尺寸 与 矩形 管道 间隙 得 到 变 
化 的 空 泡 份额 ， 而 其 余 四 种 模型 则 认为 临界 空 泡 份 
额 保持 不 变 ， 因 此 临界 空 泡 份额 的 计算 是 泡 状 流向 
单 状 流转 换 准则 的 关键 。 

表 1 模型 平均 绝对 误差 及 平均 误差 比较 (单位 : 9) 


Table 1 Mean deviation and average deviation 


tuu 


calculated for different correlations (unit: %) 

通道 尺寸 (mm) Ishii Jones-Zuber Taitel | Sadatomi 
s=0.3, w-12 24.7 29.6 50.3 53.4 
s=0. 6, w-12 66.6 70.8 81.5 83.1 
s-l. w-12 52.9 65.3 66.2 69.3 
s=1.6, w-43 10.3 13.5 61.2 66 
s=2.4, w-40 27 46.5 29.9 33.6 
s=3, w-43 43.6 18.7 50 56.4 
s=5, w-40 64.1 59.4 79.6 95 
5-10, w-43 62.4 84.9 99 100 


3. 3 泡 状 流 一 一 弹 状 流转 换 准 则 修正 

泡 状 流向 弹 状 流转 换 的 原因 ， 经 典 理论 认为 是 
气泡 间 的 相互 碰撞 、 合 并 导致 气泡 的 尺寸 增加 ， 以 
至 于 超过 通道 的 直径 而 形成 气 弹 的 结果 ， 因 此 气泡 
的 初始 尺寸 对 于 流 型 转换 的 临界 空 泡 份额 是 有 很 大 
IWD, Ishi 模型 在 窄 人 第 形 通 道中 的 适用 性 很 


结果 与 预测 结果 误差 很 小 ， 因 此 在 较 大 间隙 的 算 
通道 中 可 以 采用 公式 〈9) 计算 临界 空 泡 份额 。 


NS 


一 


液 相 折算 速度 aus) 


02F 


通 35 —f 
38 道 三 实验 结果 
通 道 三 预 测 结 果 


HUS -— (m/s) ! 
5 实验 值 与 预测 值 比较 


Fig. 5 Comparison between experimental data and new 


D 


correction 


4 结论 T6 


(OD 基于 漂移 流 思 想 ， 对 三 个 实验 段 的 分 布 
系数 进行 计算 ， 结 果 表 明 在 窗 矩 形 通道 中 可 以 选用 
Ishii 提出 的 分 布 系数 计算 式 ， 而 常规 矩形 通道 中 分 
布 系数 为 1.2。 

(2) 抱 形 通道 间隙 大 于 2 mm 时 ， 随 着 通道 水 
力 直径 的 增 大 ， 转 换 边 界 左 移 。 

(3) 在 罕 通 道 范 围 可 以 采 
流 公式 作为 流 型 转换 边界 。 


] Ishii 提出 的 漂移 


(4) 常规 矩形 通道 中 ， 转 换 准 则 中 的 临界 空 
泡 份 额 与 气泡 初始 尺寸 有 关 ， 使 用 赵 建 福 计 算 临 界 
空 泡 份额 的 计算 式 修正 Ishii 计算 空 泡 份额 的 公式 ， 
转换 边界 的 预测 与 实验 结果 符合 程度 较 好 。 
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